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Abstract 

The pairwise quantum entanglement of sites in disordered electronic one-dimensional systems 
(rings) is studied. We focus on the effect of diagonal and off diagonal disorder on the concurrence 
Cij between electrons on neighbor and non neighbor sites i,j as a function of band filling. In the 
case of diagonal disorder, increasing the degree of disorder leads to a decrease of the concurrence 
with respect to the ordered case. However, off-diagonal disorder produces a surprisingly strong 
enhancement of entanglement. This remarkable effect occurs near half filling, where the concurrence 
becomes up to 15% larger than in the ordered system. 

PACS numbers: 03.67.Mn,71.23.An,71.23.-k 
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I. INTRODUCTION 



It has been shown that entanglement is a necessary requirement for the reahzation of 
several algorithms used in quantum computation and information. For example, quantum 
teleportation, which is a technique for moving a quantum state between two points, requires 
the use of entangled states^. On the other hand, there are some evidences that many 
problems involving strong-coupled many body systems can be characterized and classified 
by using quantum entanglement. Therefore, this can lead to new and more effective methods 
for understanding the dynamical behavior of complex quantum systems^. In particular, 
it has been conjectured that entanglement can play an important role in quantum phase 
transition o^i^i^i^ . 

Several relations between entanglement and quantum phase transitions have been ob- 
tained recently. Osterloh et al., have analyzed the behavior of entanglement near the critical 
point of the spin 1/2 model XY in a transverse field^. They have found that in the region 
close to a quantum phase transition the entropy of entanglement obeys a scaling law. Al- 
though the entanglement is not an indicator of that phase transition, they showed that there 
exists an intimate connection between entanglement, scaling and universality. Moreover, Vi- 
dal et al.— , have established a precise connection between concepts of quantum information, 
condensed matter physics, and quantum field theory by showing that the behavior of critical 
entanglement in spin systems is analogous to that of entropy in conformal field theory. In 
addition, Gu et atf- have studied the properties of entanglement in the ground-state of an 
antiferromagnetic XXZ chain. They have shown that the competition between quantum 
fluctuation and ordering leads to a maximization of entanglement at the isotropic point. 

Interestingly, also in the case of electron systems a quantitative analysis of entanglement 
can be performed. Recently, Zanardi et al. studied the pairwise entanglement between 
nearest-neighbor sites in itinerant spinless-fermion systems described by a periodic tight- 
binding Hamiltoniani. They calculated the concurrence, which is a quantity which provides 
the same information as the entropy of entanglement^, and showed that the maximal 
concurrence between two neighboring sites is obtained at half band filling x = n/Na = 0.5, 
being the n the number of electrons and Na the number of atoms (sites). The concurrence 
is symmetric with respect to the half-filling point. 

The impurities or disorder effects have been seen in quantum computation as something 
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not desired that can induce errors in gating or in the evolution of a quantum state-^^. In 
addition, if the imperfection strength is increased, new phenomena related with chaotic 
behavior can occur—. On the other hand, it is well known that crystalline systems with 
impurities can undergo a metal-insulator transition as a function of disorder strengtbi^. 
Therefore, it is interesting to study the effect of disorder on the concurrence and if it is 
possible to find a relation between both. 

In this paper, we present the first study of entanglement in disordered electronic systems. 
We analyze the effect of disorder on the entanglement between sites in one- dimensional 
rings. For this purpose, we consider diagonal and off-diagonal disorder and different degrees 
of disorder. We show that whereas diagonal disorder reduces and suppresses entanglement 
for all values of the band filling x, off-diagonal disorder considerably enhances the entropy 
of entanglement for x near, but not exactly at, half-filling. 

The paper is organized as follows. In section |nj we discuss the employed model Hamil- 
tonian and how to quantify the entanglement. In Sec. IIIH the effect of the length of the ring 
on the entanglement is studied. We analyze the effect of diagonal and off-diagonal disorder 
on entanglement in Sec. IIV Al and HVBl Finally, in Sec. El we summarize our results. 

II. THEORY 

We consider one-dimensional disordered electronic systems and describe them by a 
Anderson-Hamiltonian of the foim^ 

H = Y1 ^i^i + H ^ijctcj. (1) 

For simplicity, we consider spinless electrons. In Eq. (Q) cf (q) is the usual creation (anni- 
hilation) operator of a spinless electron at site z, whereas — Cj IS the number operator, 
and tij is the hopping integral between nearest neighbor (NN) sites i and j. Si is the on-site 
energy for atom i. Diagonal disorder is introduced by assuming that the values of Si are 
randomly distributed in the energy interval [0, W] while the NN hopping integrals remain 
unchanged with respect to the perfectly ordered case, i.e., Uj = —1. Analogously, the effect 
of off-diagonal disorder (random hopping) is modeled by considering the hopping integrals 
tij to acquire random values in the interval [0, W] while all on-site energies are equal to zero, 
like in the ordered lattice. 
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To calculate the entanglement of formation between pairs of sites we resort to the concept 
of concurrence^!^ and proceed in a similar way as Zanardi and coworkers did for the case of 
translationally invariant chains^. 

Basically, one can consider the ring with A^^ sites and n electrons as being a bipartite 
system composed of two subsystems, A and 5, where A refers to the pair of sites i and 
j and B to the rest of the ring (see Fig. 1). More precisely, A contains all 4 possible 
electronic states in the pair of sites i and j, i.e., \1)a = \0)a, |2)a = cf\0)A, |3)a = c^|0)yi, 
and \4:)a = cfcj'\0)A- Analogously, B contains the states of the rest of the ring. Thus, a 
particular state I'^ab) of the entire ring can be represented as a direct product of the form 

\^ab) = E ® li^Bu), (2) 

where \'^pA^l) and \ipBu) are basis states of the subsystems A and B. Thus, the density matrix 
of the whole system reads 

p= (3) 

The reduced density matrix pA for the subsystem A is obtained by performing a trace on 
subsystem B in order to integrate out its degrees of freedom. Therefore 

PA = TrB{\^AB){^AB\}. (4) 

The entropy of entanglement of a bipartite system is calculated as^ 

EH'^ab)) = -Tta {pA logpA} = -Ttb {pB logpfi} (5) 

and quantifies the entanglement of formation. It has been shown by Wootters^ that the 
entropy of entanglement can be also characterized by the concurrence C, since E^Y^ ab)) is 
a monotonous function of C, which is defined by 

C = max (Ai — A2 — A3 — A4, 0) . (6) 

The Afc's are, in decreasing order from = 1 to A; = 4, the eigenvalues of the matrix 
R = VPa Pa, where pa = {cry ® cXy) p\ {ay ® ay), where p\ is the complex conjugated of the 
matrix pA and ay the "y" Pauli-matrix. 

Thus, in order to calculate C for disordered ID systems described by the Hamiltonian of 
Eq. (dl), we should first determine pa by means of Eqs. (j2I) and (jH). As pointed out in Ref. 
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the number of electrons n is conserved due to the fact that [H , n] =0. As a consequence, 
only the elements p^(ll), pa(22), p^(33), pa(44), pa(23) and pa(32) are nonzero. All other 
elements vanish. This can be seen by looking at the expansion When constructing 
\^ab){^ab\ and performing the trace over B one obtains, for all elements of pA others 
than the ones mentioned above, scalar products of the type (ipBuli^Bfj.) , between states of B 
having different number of particles and, therefore, being orthogonal. One can take profit 
of this property of the reduced density matrix and obtain its nonzero elements as quantum 
expectation values of relevant physical quantities. For this purpose we use the whole wave 
function |\E')a_b of the system and write: 

Pa(11) = {^AB\{l-n^){l-n,)\^AB) (7) 

PAiU) = {^ABlflifljl'^AB) (8) 

PAi22) = {^AB\n^il-n,)\'^AB) (9) 
Pa(33) = - ni)nj\^AB) (10) 

p^(23) = p:^(32) = (^abIc+c.I^ab). (11) 

In Eq. ((7j), the presence of the operator (1 — hi){l — fij) automatically projects out all states 
of the whole system not having electrons on site i and electrons on site j. Similarly, the 
operator -hifij in Eq. (jH)) ensures that only states having 1 electron on i and 1 electron on j 
will contribute to the element p^(44). In Eqs. Q and (fTI]|) the corresponding projectors to 
obtain pyi(22) and pyi(33) are described. Finally, pa(23) [pa(32)] contains only contributions 
from states having 1 electron on site i (j) and no electrons on site j (i). 

Thus, with the help of the above equations, one can write reduced density matrix pA in 
the following iorm-^ 



PA 



fv 0\ 

w z 
z u 



(12) 



\0 yj 
where the nonvanishing matrix elements are given by 



V = 1 - {fii) - ijij) + {fiifij), y = {fiihj) 

z = {c^Ci), w = {fii) -y, u= {fij) - y. (13) 

Eqs. (fT^ show how the nonzero elements of pA, can be calculated as average quantities of 
the complete ground-state wave function. 



Diagonalization of the matrix of Eq. ()12p is trivial. Using Eq. © one obtains an analytic 
expression for the concurrence between sites i and j of the forn>^ 

Cij = 2max{0,\z\-^}. (14) 

Eq. (fT^ does not contain w and u. Note that i and j do not need be nearest neighbors. In 
the next section we present calculations of Ci^i+e with i = 1,3. 

For every set of parameters, the Hamiltonian can be diagonalized numerically and the 
concurrence C can be obtained using Eqs. (fT^ and (fT^ . The average value of the con- 
currence can be then calculated by using a large number of replicas in order to minimize 
numerical fluctuations. 

To calculate the average concurrence of the tight-binding rings with diagonal and off- 
diagonal disorder, we have implemented the following procedure: (i) we have generated a 
set of random parameters (on-site energies or hopping elements) for each of the replicas in 
order to model the disorder of the finite size rings, (ii) we have numerically diagonalized the 
Hamiltonian, given by Eq. Q for obtaining its single-particle eigenvectors, [iii) with the help 
of these eigenvectors, we have calculated the values v, y and z necessary for the calculation 
of the concurrence and, finally, (iv) we have obtained the average concurrence Cij {Cu, C13, 
6*14) as 

1 ^« 

^^■1+^ ^ N — '■ — iV~ ^ ^ (15) 

replicas a replicas k 

We have performed these calculations for different ring sizes. We found that for rings 
consisting of Na — 200 sites results achieve convergence. Therefore, in this paper all the 
results for Cij correspond to A^^^ = 200 sites, unless we explicitly specify the ring size. On the 
other hand, the number of replicas we used was 10000. In order to understand the influence 
of disorder on the concurrence it is necessary to study first the behavior of the concurrence 
in ordered systems. 

III. ORDERED RINGS: SIZE EFFECTS 

The concurrence between C12 as a function of n for a perfectly ordered, periodic ring can 
be easily computed, using the theory presented in the previous section, as a limiting case of 
Cij for Ei = OWi and tij = —1 ). Obviously, the sums of Eq. (fT3j) are not necessary in this 
case, since they cancel with the factors in the denominator. We have computed C12 vs n for 
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TABLE I: Calculation of the elements of the reduced density matrix and the concurrence 
according to equations to , for ring consisting of 4 atoms and n electrons. 



n 


1 


2 


3 


4 


\ J- / 


1/4 


3/4 


3/4 


1 




1/4 


1/4 


3/4 


1 


Z 


1/4 


1/4 


1/4 





(nin2) 





1/8 


1/2 


1 


V 


1/2 


1/8 








y 





1/8 


1/2 


1 


Cl2 


1/2 


1/4 


1/2 






different sizes. Results are shown in Fig. |21 As an example, we show below the analytical 
calculation for = 4. 

The 4x4 Hamiltonian has the eigenvalues ei = —2, £2 = £3 = 0, 64 = 2, and the 
eigenvector matrix reads 





1 





*\ 


1 
2 





1 


1 

2 


1 
2 


1 

V2 





1 

2 







1 

V2 


4/ 



(16) 



With this information at hand it is possible to calculate the elements of the reduced 
density matrix and therefore the concurrence, as summarized in Table I. 

Clearly, the concurrence is symmetric with respect to n = 2 (half filling). The concurrence 
for n = and n = 4 is equal to zero. We observe in Fig. |2ta) that for smaller ring {Na = 4, 
8 and 12 sites) it is very difficult to perceive which is the band filling x where Cu will have 
a maximum. In addition, the concurrence as a function of n does not have a monotonous 
behavior. These size effects can be understood with the help of the correlation functions 
z = {C1C2), y = {nin2), and the band filling x. It is easy to show that z = {C1C2) is related 
to the ground-state energy Egs, i.e. z = {C1C2) = —Egs/^Na, while that v = (x — 1)^ — z^ 
and y = x'^ — z"^ (see Ref. 7). From these expressions, we obtain for systems with band filling 
X = 1/Na (only one electron in the lattice) that C12 = 2z = —Eg^/Na where y = (^1^2) = 0. 
This explains the value of Cu found for Na = 4: with one electron {x = 0.25). In general, 
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the concurrence of periodic rings can be written as follows^ 



Ci2 = 2max|o, \Ej2Na\ - ([(x - 1)' - {E^,/2Naf] [x^ - {E^s/2Na) 



1/2 



(17) 



When we begin to increase the ring size, we notice that the finite size effect begins to 
dissapear [see Fig. |2fb)]. Moreover, C12 begins to have a monotonous behavior and we find 
that the maximum value of NN concurrence occurs at half-band filling. For Na = 64 sites 
the thermodynamic limit has been reached and the values of C12 are almost converged. 



IV. DISORDERED RINGS 
A. Diagonal disorder 

The electronic properties of strongly disorder systems have been investigated by a variety 
of analytical and numerical methods. In this paper, we study this problem from a very 
different perspective, the perspective of the concurrence. In Fig. El we show the NN concur- 
rence C12 for a ring with Na = 200 sites where we have introduced only diagonal disorder. 
We observe that the effect of this diagonal disorder is to decrease C12 in a monotonous way. 
In Fig int^a) C12 is presented as a function of band filling for some representative values of 
the disorder strength W. This decrease of the concurrence is related to the localization of 
the wave function. When disorder is introduced, for example on site i, the wave functions 
are localized around this site, therefore the correlation function z = (cfcj) between two NN 
sites is decreased. On the other hand, Osterloh et al.— , in their study of the critical point in 
the XY model in a transversal magnetic field, showed that the behavior of the concurrence 
does not have a divergence in the transition point. However, they showed that the derivative 
of C12 shows this divergence and scaling properties. Therefore, it may be of interest to study 
the derivative of C12 as a function of the disorder strength. Unfortunately, it is known that 
numerical calculation of derivatives in disorder systems as a function of W displays errors 
associated with the numerical divisioni^. In order to visualize the behavior of the derivative 
of C12 as a function of the strength of disorder, in Fig Ol^b) we plot C12 as a function of the 
disorder strength W for some representative values of the band filling x. From the figure, 
we can deduce that the derivative of C12 as a function of disorder strength is almost zero 
in the range of G [0, 2.5]. After that, the derivative of C12 as a function of W is almost 
constant with a slight negative slope. 
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A similar analysis can be realized on the concurrence as a function of the distance between 
sites. In particular, we have focused on next nearest-neighbor concurrence C13 and third 
nearest-neighbor concurrence C14. The results are shown as a function of the band filling x 
for some representative values disorder strength W [Figs. El^a) and Eta)] and as a function 
of the disorder strength W for some values of the band filling x [Figs, ^b) andEfb)]. We 
observe that the values of the concurrence for all band fillings and disorder strengths is 
much smaller than Cu- This implies that in systems without disorder and systems with 
diagonal disorder the concurrence is highly local and strongly decreases as a function of 
the distance between sites. Notice that the maximum value of the concurrence occurs at 
W = and begins to move towards band fillings smaller than the half band 1/2 when we 
increase the distance between sites and disappears in some range of band filling. This result 
can be explained by using the fact that the correlation function z = (cfcj) depends on the 
overlap between the wave function, which decreases when we increase the distance between 
sites, whereas v and y fundamentally depend on the band filling. Therefore, it will have 
a threshold band filling where ^/vy > \z\ and, consequently, Cij = 0. On the other hand, 
as we begin to increase the disorder strength, the concurrence is almost insensitive to the 
system band fillings. These behaviors can be observed in figures 4(b) and 5(b), where we 
show the concurrence for systems with different band fillings and W > 10 are joined in the 
same curve. This result is related to the fact that due to the strong disorder z, v are y small 
for all band fillings. 

B. Off diagonal disorder 

Although for ID systems with diagonal or off-diagonal disorder the wave function are 
localized for all disorder strengths, it is known that the physical properties are strongly 
dependent on the type of disorder. Systems with off-diagonal disorder show an anomalous 
behavior in the density of states and in the localization lengtbi^*ii, where it has been found 
that E = state is localized but has an infinite localization length. In this subsection, we 
show how this property affects the concurrence. 

In Fig. ini we present results for Cu for a ring with Na = 200 sites and with off-diagonal 
disorder ^nn £ [0, W]. In Fig. Efa) ^712 is showed as a function of band filling x for some 
representative disorder strengths W whereas in Fig. IH^b) ^712 is presented as a function 
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of W for some representative values of x. From the figures, we observe that, contrary to 
the diagonal disorder, C12 value can increase for some band fillings in comparison with the 
non disordered case {W = 0) when the disorder strength is increased. In particular, we 
observe that this occurs for band fillings larger than x > 0.25 and for off-diagonal disorder 
strength W > 4. Notice that the maximum value of the concurrence is found for band filling 
values in the range 0.4 < x < 0.45 and for W > 20. Moreover, the overall form of the 
concurrence curves does not change anymore for W > 20. It is surprising that in the case of 
system with off-diagonal disorder, the NN concurrence increases in comparison with respect 
to the non disordered case. This behavior can be related to the anomaly in the density 
of states found in these kind of systems. To inquire about it, we have calculated the NN 
concurrence as a function of the band filling for a small system with A^^ = 16 sites with 
only one different hopping elements with respect to the others. These results are shown in 
Fig. [71 where we present the concurrence between the sites with the different hopping integral 
(ti2 = — 2,tNN = —1, otherwise), Cu, and the NN concurrence between sites different from 
1 and 2. In the last case, we only show the concurrence between sites close to the impurity. 
From the figure we observe that C12 has much larger values compared to the result show 
in the previous subsection. This concurrence is almost independent of the band filling and 
oscillates around 0.72. These high values of C12 are related to the fact that the impurity 
localizes the eigenfunction with the lowest eigenvalue Eq between the sites where the 
impurity was placed (ti2 = —2), whereas the other eigenfunctions are almost delocalized. 
Thus, we have that z = {ipolcf C2\'4'o) has very large contribution to the concurrence for all 
band fillings. On the other hand, we know that in the case of one electron, vy = 0. Therefore 
from the figure, we can see that the quasidelocalized eigenfunctions contribute very little 
to the NN concurrence when we increase the band filling. In the case of C23, we observe 
that the localization effect turns C23 into a decreasing function. It is necessary to remark 
that the localization effect due to one impurity only strongly modifies the NN concurrence 
between sites very near to this, where the eigenfunction \tpo) is localized, whereas the NN 
concurrence of remote sites from the impurity show a behavior more similar to the case of 
a periodic ring. Therefore, when the hopping integrals are randomly chosen in the range 
i^NN £ [0, W^], the number of hopping integrals having values larger than W/2 is Na/2. 
This implies that in general there are Na/2 eigenfunctions that are localized on the bonds 
affected by these hopping integrals. These Na/2 bonds will contribute with very high values 
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to Cnn, as has been shown in the case of one impurity, whereas the other Na/'i bonds will 
contribute with much smaller values. As we increase W , the localization of the eigenfunction 
will be increased around these bonds and their C^^^ values will be close to 1 whereas the 
Cnn values of the other bonds will be close to zero. These bond-localizations explain the 
increase in the NN concurrences of systems with off-diagonal disorder in comparison with 
the non disordered case. In addition, this clearly shows the difference between both types 
of disorder, in diagonal disorder the one particle eigenf unctions are localized around the 
sites, which decreases the total value of Cnn with increasing W . On the other hand, in the 
off-diagonal disorder case the eigenfunctions are localized between the nearest neighbors, 
which increases Cnn with increasing W . 

In order to study the way the off-diagonal disorder affects the concurrence as a function 
of the distance between sites in this kind of systems, we have calculated the concurrence as 
a function of the distance. The results are shown in Figs.|Ht'nd|Hl In these figures we observe 
that, as the off-diagonal strength W is increased, the maximum in the concurrence begins 
to decrease, contrary to the case of the NN concurrence. This confirms that the effect of the 
increase in the maxima of the NN concurrence is due to the strong localization on the bonds. 
Contrary to the case of diagonal disorder, the decrease of the maximum values of Ci^ and 
6*14 for large W is relatively small with respect to the non disordered case. We again notice 
that the concurrence decreases as a function of the distance between sites, which shows that 
the concurrence is a strongly local property. Moreover in these cases, the curves of Ci^ and 
6*14 depend on the band filling. This behavior can be seen in Figures 8(b) and 9(b), where 
we show that the curves for different band fillings are completely different. Finally, for a 
given band filling x and for W > IQ, the curves are almost insensitives to the strength of 
the disorder. 

V. CONCLUSIONS 

In this work, we have studied the effect of diagonal and off diagonal disorder on ID 
electronic systems described by the Anderson Hamiltonian. We have found that, in the 
case of diagonal disorder, as we increase the disorder strength the concurrence between 
sites decreases in a monotonous way. In addition, we found that the concurrence in these 
systems is highly local and decays for increasing distances between sites. On the other 
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hand, in off diagonal disorder case we observe that the increase in the disorder strength 
results in an increase of the NN concurrence in comparison to finite rings in the band filling 
interval 0.25 < x < 0.5. This behavior of the concurrence for both kinds of disorder can be 
understood with the help of the one particle wavef unctions. In the diagonal disorder case, 
the one particle wave functions begin to localize on the sites as we increase the disorder 
strength. Therefore, this localization decreases the NN concurrence. On the other hand, 
the increase of the off-diagonal disorder strength localizes the one particle wavefunction on 
bonds between NN sites, which increases the NN concurrence in a certain interval of the 
band filling. This result could be useful for the study of entanglement in 3D amorphous 
solids. Although the metal-insulator transition occurs only in 3D systemsi^, it is known 
that the off-diagonal disorder in ID systems produces anomalous propertie o^^i^^ that should 
have some effect of the concurrence. 

Investigations of the concurrence as an indicator of the insulator-metal phase transition 
in 3D disordered systems are in progress. 
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FIG. 3: (Color online) Nearest neighbor concurrence average C12 of a ring with Na = 200 sites . In 
(a) we present C12 as a function of band filling x for several representative values of the diagonal 
disorder strength W. In (b) C12 is plotted as a function of disorder strength W for some band 
filling X. -j^g 
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FIG. 4: (Color online) Next nearest-neighbor average concurrence C13 of a ring with Na = 200 
sites as a function of band filling x for some representative values of diagonal disorder strength W. 
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FIG. 5: (Color online) Third neighbor average concurrence (7x4 of a ring with Na = 200 sites as a 
function of band filling x for some representative values of diagonal disorder strength W 
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FIG. 6: (Color online) Nearest neighbor average concurrence C12 of a ring with Na = 200 sites as 
a function of band filling n for several representative values of the off diagonal disorder strength 
W. 
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FIG. 7: Nearest-neighbor concurrence Cu for a ring with Na = 16 sites and one impurity. The 
impurity is locahzed between the sites 1 and 2 {tu = —2). The hopping integral between the other 
nearest neighbor sites is iNN = — 1.0 
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FIG. 8: (Color online) Next nearest neighbor average concurrence C13 of a ring with Na = 200 
sites as a function of band filling n for several representatives values of the off diagonal disorder 
strength W. 
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FIG. 9: (Color online) Third neighbor average concurrence C14 of a ring with Na = 200 sites as a 
function of band filling n for several representatives values of the off diagonal disorder strength W . 
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